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SECTION  I 


INTRODUCTION 


One  of  the  unique  characteristics  of  laminated  composite  materials 

is  the  delamination  mode  of  failure  that  occurs  at  free  edges.  This  behavior 

has  been  investigated  experimentally  and  by  mathematical  stress  analysis 

models  for  laminates  with  straight  free  edges.  Experimental  studies  by  Pipes, 

[11* 

Kaminski,  and  Pagano  have  shown  that  the  laminate  stacking  sequence  can 

effect  the  static  strength  of  laminates.  Similar  effects  of  stacking  sequence 

[2] 

were  found  by  Foye  and  Baker  for  fatigue  loading  conditions.  Mathematical 

stress  analyses  of  laminated  plates  with  straight  free  edges  have  shown  that 

changes  in  stacking  sequence  can  alter  the  behavior  of  the  interlaminar  stresses. 

In  particular,  investigators ^»4,5]  have  found  that  the  sign  of  the  transverse 

normal  stress  o changes  from  tension  to  compression  as  the  stacking  sequence 
Z r 61 

is  changed.  Pagano  and  Rybicki1  fotind  that  a change  in  stacking  sequence 

of  the  type  (0°/Matrix),  to  (Matrix/0o)  reversed  the  sign  and  changed  the 

® f 7 8 1 

magnitude  of  the  interlaminar  stress,  a^.  Furthermore,  Pagano  and  Pipes1  * 

have  shown  that  high  tensile  stresses  are  associated  with  the  decreased  lamin- 
ate strengths  reported  in  References  [1]  and  [2].  These  observations  point 

to  the  importance  of  understanding  the  stress  behavior  of  laminates  near  a 

[9] 

free  edge.  Through  such  an  understanding,  Lackman  and  Pagano  have  success- 
fully reduced  interlaminar  stresses  at  straight  free  edges  to  suppress  the 
delamination  mode  of  failure. 

While  the  straight  free  edge  can  be  studied  as  a problem  with  two- 
dimensional  variations  in  stresses  and  displacements,  the  curved  free  edge 
is  inherently  a three-dimensional  problem.  Thus,  the  circular  hole  problem 
is  more  difficult  to  treat  in  terms  of  a mathematical-stress  analysis  than 
the  straight  free-edge  problem.  Dana  and  Barker^1^  have  investigated  the 
three-dimensional  stress  analysis  of  four-ply  (0/90)  and  (45/-45)  laminates 
with  circular  holes  and  found  that  the  sign  and  magnitude  of  the  o^  distribution 
around  the  hole  can  be  changed  by  altering  the  stacking  sequence.  Rybicki  and 
Hopper^1^  have  developed  an  analysis  for  the  interlaminar  stress  distributions 


* Numbers  in  brackets  denote  references  given  at  the  end  of  the  report. 
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for  a six-ply  symmetric  laminate  containing  a circular  hole.  Experimental 

f 12] 

investigations  by  Daniel,  Rowlands,  and  Whiteside  on  the  effects  of  stack- 
ing sequence  on  the  strength  of  laminated  plates  with  holes  showed  that  varia- 
tions in  stacking  sequence  effected  the  laminate  strength  and  failure  mode. 

Other  studies  for  example  by  Greszczuk^^ , Waddoups,  et  al.^^, 
Waszczuk  and  Cruse  , and  Whitney  and  Nuismer^^  have  predicted  failure 
of  a composite  plate  containing  a hole  and  made  comparisons  with  experimental 
data.  References  [13]  and  [14]  are  concerned  with  predicting  failure  for 
damage  due  to  inplane  stresses.  The  stress  analyses  for  these  studies  did 
not  include  interlaminar  stresses  and  thus  are  not  directly  applicable  to  the 
delamination  problem. 

There  are  several  basic  studies  that  can  provide  insight  into  under- 
standing laminate  behavior  near  a curved-free  edge.  Five  such  studies  are 
considered  here.  The  first  study  concerns  the  effect  of  stacking  sequence 
on  the  transverse  stresses.  For  the  problem  of  the  straight  free  edge. 
References  [3]  through  [9]  point  out  that  changing  the  stacking  sequence  can 
change  both  the  sign  and  magnitude  of  the  interlaminar  stress  distribution. 
However,  in  the  case  of  a laminate  containing  a circular  hole,  the  effect  of 
changing  the  stacking  sequence  is  not  understood  as  well.  The  second  study 
considers  the  influence  of  lay-up  angle  on  the  interlaminar  stress  around 
circular  holes.  The  third  study  pertains  to  the  state  of  residual  stresses 
that  exist  in  the  laminate  due  to  the  temperature  changes  in  the  fabrication 
process.  In  the  fourth  study,  the  lnplane  stress  distributions  around  the 
circular  hole  obtained  from  a three-dimensional  model  are  compared  with  two- 
dimensional  model  results  obtained  from  a laminated  plate  theory  model.  In  the 
fifth  study,  the  effect  of  laminate  thickness  on  transverse  stresses  for  a 
plate  with  a given  hole  size  is  examined  for  a selected  laminate.  This  prob- 
lem is  of  primary  importance  in  the  mathematical  modeling  of  laminated  plates 
because  a thin  plate  with  a large  hole  can  present  high  element  aspect  ratios 
or  a large  number  of  elements  in  a finite  element  analysis.  In  addition,  a 
comparison  of  predicted  and  measured  strains  around  a circular  hole  and  an 
approximate  method  for  predicting  the  sign  of  the  interlaminar  stress,  a z> 
are  presented.  These  topics  are  an  important  step  toward  obtaining  a better 
understanding  of  the  behavior  of  laminates  containing  a hole  and  are  pertinent 
to  designing  laminated  structural  components  as  well  as  to  designing  test 
specimens  and  Interpreting  laboratory  test  results. 
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SECTION  II 


METHOD  OF  STRESS  ANALYSIS 


The  method-of-stress  analysis  is  based  on  a compatible  finite  element 
representation  for  the  deformations  of  a three-dimensional  solid.  The  computer 
program  is  SAP  IV.  Details  of  this  program  and  the  stress  analysis  method  are 
given  in  Reference  [17].  A summary  of  the  analysis  is  given  here  as  back- 
ground information. 

In  the  finite  element  stress  analysis  plys  are  modeled  as  three- 
dimensional  homogeneous  orthotropic  materials  with  linear  stress-strain  behavior. 
Symmetric  laminates  are  considered  with  applied  stress  loading  conditions  dis- 
tant to  the  hole.  The  analysis  is  based  on  a three-dimensional  finite  element 
representation  for  the  displacements  within  each  element.  Compatibility  con- 
ditions within  each  element  and  between  contiguous  elements  are  satisfied  ex- 
actly. Equilibrium  is  satisfied  approximately.  A twenty-node  element  with 
corner  and  midside  nodes  was  used. 

Material  properties  are  based  on  effective  modulus  representations 
for  the  plys.  This  representation  does  not  recognize  fiber  size  and  geometry. 

However,  effects  of  such  characteristics  can  be  evaluated  through  the  approach 

r ^g  I 

described  by  Rybicki  and  Pagano1  . 


3 


jL 


SECTION  III 


NUMERICAL  RESULTS 


This  section  contains  numerical  results  for  cases  selected  to  address 
the  question  concerning  effects  of  stacking  sequence,  lay-up  angle,  thermally 
induced  residual  stresses,  and  laminate  thickness  on  the  stress  distribution 
around  a hole  in  a laminated  plate.  In  addition,  some  comparisons  of  three- 
dimensional  and  two-dimensional  laminated  plate  theory  solutions  are  given. 

The  numerical  results  presented  in  this  section  are  based  on  seven 
laminates  with  different  lay-ups.  Each  laminate  was  subjected  to  a mechanical 
stress  loading  and  three  of  the  laminates  were  selected  for  a thermal  loading 
condition.  In  addition,  each  of  these  ten  cases  were  run  for  a different 
stacking  sequence.  Thus,  a total  of  twenty  analyses  were  performed.  Because 
of  symmetry,  only  one  quarter  of  the  top  half  of  each  laminate  was  modeled. 

In  each  case,  a thirty-six  element  grid  with  three  elements  through  the  thick- 
ness and  twelve  in  the  laminate  plane,  as  shown  in  Figure  1,  was  used  to 
model  one  eighth  of  the  laminate. 

1.  Effect  of  Stacking  Sequence 

The  effect  of  stacking  sequence  has  been  an  important  topic  of  study 
for  free-edge  problems.  As  stated  previously,  the  effects  of  stacking  sequence 
on  a circular  free  edge  are  not  as  well  understood  as  for  the  case  of  the 
straight  free  edge.  To  provide  more  insight  into  this  problem,  seven  laminates 
with  a mechanical  stress  loading  condition  and  three  laminates  with  a thermal 
stress  loading  condition  were  considered.  Results  were  also  obtained  for  the 
associated  ten  cases  with  different  stacking  sequences. 

The  lay-up  for  the  first  laminate  was  (0/90)g.  Ply  properties  for 
the  0°  ply  are  given  in  Table  1.  The  applied  stress  boundary  conditions  were 
selected  to  produce  a unit  average  stress  parallel  to  the  0°  direction  or  in 
the  y direction  illustrated  in  Figure  1.  The  applied  stresses  are  given  in 
Table  2.  The  convention  for  lay-up  angles  is  that  angles  are  measured  clock- 
wise from  the  load  direction.  A laminate  thickness  to  hole  diameter  ratio  of 
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TABLE  1. 

MATERIAL  PROPERTIES 
LAMINA  IN  THE  (0/90) 

OF  THE  0" 

LAMINATE 

s 

• 

Extensional 

Modulus 

(psi) 

Poisson's 

Ratio 

Shear 

Modulus 

(psi) 

E1 

= 30.0  x 10^  psi 

v12  = .336 

G12 

- 1.0  x 

io6 

E2 

= 3.0  x 10^  psi 

v13  " *336 

G13 

= 1.0  x 

106 

E3 

= 3.0  x 10^  psi 

v23  = ,336 

G23 

- 1.0  x 

h- 

O 

O' 

TABLE  2.  STRESS 
(0/90)s 

BOUNDARY  CONDITIONS  FOR 
AND  (90/0)  LAMINATES. 

Material 

°x(psi) 

Oy(psi) 

T (psi) 
xy 

0° 

.05 

1.821 

0. 

90° 

0 

1 

.179 

0. 

d 


* 


.8  was  selected.  The  results  for  the  c ^ distribution  around  the  circular  hole 

at  the  laminate  midplane  are  shown  in  figure  2.  Two  Important  results  are  shown 

in  this  figure.  First,  the  stress  distribution  changes  sign  and  magnitude  as 

the  stacking  sequence  is  changed.  Second,  the  sign  of  at  8 * 90  for  these 

(90/0)  and  (0/90)  laminates  can  be  predicted  from  the  sign  of  the  moment  away 
s s 

from  the  circular  hole  along  0 m 90,  by  using  equilibrium  considerations  on 

the  remote  applied  stresses,  but  this  is  not  true  at  6 = 0°  where  the  sign  of 

the  moment  due  to  remote  stresses  would  always  indicate  a tensile  stress  for  a^. 

Six  other  laminates  were  considered.  These  were  of  the  type  (+0/+0/O2)8 

and  (+0/+0/9O„)  where  0 is  30",  45°,  and  60°.  Stacking  sequences  were  changed 
— Zb 

by  reversing  the  (+0/+0)  plys  and  the  (0^)  or  (90j)  plys.  The  (0^)  and  (9O2) 

plys  were  modeled  as  one  material  and  the  (+0/+0)  plys  were  modeled  as  a single 

material  with  effective  modulus  properties.  Table  A-l  of  Appendix  A contains 

the  properties  obtained  from  Reference  [19]  and  used  in  this  part  of  the  study. 

A unit  stress  in  the  0°  direction  was  applied  to  each  laminate.  The  stress 

boundary  conditions  for  each  laminate  are  given  in  Table  A-2  of  Appendix  A.  In 

addition  to  the  unit  stress  loading,  the  effect  of  stacking  sequence  for  a thermal 

loading  of  -200  F was  also  considered.  Coefficients  of  linear  thermal  expansion 

shown  in  Table  A-l  were  used  in  this  calculation.  For  the  thermal  loadings,  only 

the  last  three  laminates  designated  in  Table  A-2  were  considered  because  in  terms 

of  the  models,  the  last  three  laminates  are  equivalent  to  the  first  three  rotated 

by  90°.  For  both  mechanical  and  thermal  loading  conditions,  a hole  diameter  of 

1 inch  and  laminate  thickness  of  .6  inch  was  used. 

Tables  A-3  through  A-8  of  Appendix  A contain  stresses  at  the  laminate 

midplanes  for  a unit  stress  applied  parallel  to  the  0°  direction.  Reference 

will  be  made  to  these  tables  to  describe  some  of  the  characteristics  of  the 

stress  distribution  for  the  interlaminar  stress  denoted  by  o . Comparison  of 

z 

Tables  A-3  and  A-6  and  Tables  A-4  and  A-7  show  that  for  the  (+0/+0/O2)8  to 

(O./+0/+0)  changes  in  stacking  sequence,  the  sign  and  magnitude  of  the  a stress 

Z ~~  B * 

distribution  at  the  laminate  midplane  is  changed.  It  is  noted  that  the  o z stresses 

for  the  (O./+0/+0)  laminates  are  generally  lower  than  the  values  for  the 

2 - s 

(+0/+0/O2)s  laminates. 

The  (+30/+  30/902)g  and  (+45/+  45/902)8  laminates  display  a o8  dis- 
tribution that  is  totally  tension  around  the  edge  of  the  hole  as  seen  from  the 
results  in  Tables  A-4  and  A-5.  The  effect  of  changing  stacking  sequence  is  to  change 
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this  distribution  to  totally  compression  as  seen  in  Tables  A-7  and  A-8.  The 

c>z  distribution  for  the  (+6O/+6O/9O2) g laminate  (Table  A-5)  has  a sign  change 

from  compression  to  tension  in  going  from  0®  to  90°.  The  associated  (902/+60/+60)s 

laminate  (Table  A-8),  however,  does  not  show  this  sign  change.  The  results 

show  1^-ger  values  for  in  the  (+e/+8/902)g  laminates  than  in  the  (902/+9/+8)g 

laminates. 

The  effects  of  changes  in  stacking  sequence  for  thermal  loading  is 

shown  by  the  results  in  Tables  A-9,  A-10,  and  A-ll.  Once  again,  a change  in 

both  sign  and  magnitude  are  found  for  changing  stacking  sequence.  Several 

interesting  features  of  these  results  are  noteworthy.  The  (+30/+30/90 2>  and 

the  (902/+30/+30)s  laminates  shew  a change  in  sign  of  o in  going  from  0=0° 

to  90°.  The  (+45/+45/90,)  , (90./+45/+45)  , (+60/+60/90J  , and  (90./+60/+60) 

2 s Z — s — 2 S l — s 

laminates  do  not  show  this  sign  change.  Generally,  higher  values  of  o^  occur 
for  the  (9O2)  ply  in  the  center.  The  two  laminates  with  (+60/+60)  show 
relatively  flat  distribution  for  o z as  a function  of  angle.  The  (902/+45/+45)s 
shows  the  highest  value  of  tensile  a . Finally,  it  is  noted  that  the  peak 
values  of  o z do  not  always  occur  at  0 = 0°  or  at  0 = 90°. 

2.  Effect  ci  Lay-Up  Angle 

Results  showing  the  effect  of  lay-up  angle  for  (O2/+0/+0)g,  (9O2/+0/+6) s , 

(+0/+9/O,)  , and  (+8/+0/9O-)  laminates  are  shown  in  Tables  A-3  through  A-ll.  One 
— 2 s — Zs 

obvious  effect  of  increasing  the  lay-up  angle  is  the  increase  in  the  peak  tangential 

stress  (at  0 - 90°)  in  the  center  plys  for  the  case  where  the  (9O2)  or  (O2) 

plys  are  in  the  center  of  the  laminate.  Peak  compressive  and  tensile  values 

of  for  the  (+Q/+6/0^)  laminates  are  relatively  close  together.  Peak  values 

of  tensile  stress  are  highest  for  the  (+45/445/902^  and  (+6O/+6O/9O2) g 

laminates.  For  the  laminates  with  the  (O2)  plys  on  the  outside,  the  peak  values 

of  oz  generally  Increase  for  increasing  values  of  0.  For  the  thermal  loadings 

considered,  the  effect  of  changing  lay-up  angle  is  to  change  the  character 

of  the  o distribution  around  the  hole.  This  is  seen  in  Tables  A-9  through 
z 

A-ll  for  the  (+0/+0/9O.)  laminates.  Of  the  three  (9O_/+0/+0)  laminates 

'—Zb  2—s 

considered,  the  maximum  value  of  tensile  a occurs  for  the  case  0 = 45°.  The 
’ z 

effects  of  lay-up  angle  also  change  0.  in  magnitude  and  sign  as  seen  for  the 

0 

thermal  loading  conditions  shown  in  Tables  A-9  through  A-ll. 


9 


3.  Residual  Stresses  Due  to  Fabrication 


Graphite/epoxy  laminates  are  cured  at  an  elevated  temperature. 

Cooling  to  room  temperature  can  cause  residual  stresses  In  the  laminate  due 
to  the  mismatch  in  the  coefficients  of  thermal  expansion.  To  obtain  an  estimate 
of  the  magnitude  and  distribution  of  these  stresses,  three  laminates  were  modeled 
for  a temperature  loading  condition  of  -200  F.  The  material  properties  are 
shown  in  Table  A-l.  The  boundary  conditions  were  stress  free.  Tables  A-9 
through  A-ll  show  the  stress  results  at  the  laminate  midplane  for  (+0/+9/9O2)g 
and  (9O./+0/+0)  laminates.  The  highest  tensile  value  for  is  9570  psi  and 

A *“  S 0 

occurs  in  the  (9O2)  plys  of  the  (j^0/+60/902)g  laminate.  The  highest  compressive 

sive  value  of  oQ  is  11,560  psi  and  occurs  in  the  (90_)  plys  of  the  (+30/+30/90.) 

v 2 — 2 s 

laminate.  Tensile  values  of  from  800  psi  to  1815  psi  are  predicted  for  the 
six  cases  with  the  maximum  value  occurring  for  the  (+30/430/902^  laminate. 

The  magnitudes  of  the  predicted  transverse  stresses  are  approaching  values 
large  enough  so  that  they  cannot  be  neglected  compared  to  the  transverse  strength 
of  the  laminate. 

4.  Comparison  of  LPT  and  Three-Dimensional  Stress  Distributions 

One  stress  analysis  approach  that  is  commonly  used  for  fiber  rein- 
forced plates  is  based  on  the  laminated  plate  theory  (LPT).  A description  of 
the  method  can  be  found  in  many  references,  for  example,  see  Reference  [18). 

In  the  LPT  models,  the  out-of-plane  stresses  are  zero.  LPT  analyses  are  used 
in  the  design  of  composite  structures  to  predict  their  survivability.  Thus, 
it  is  of  interest  to  know  how  the  LPT  stresses  compare  with  the  three-dimensional 
stress  distributions.  In  the  far-field  regions,  away  from  edges  and  attachments, 
the  LPT  solution  is  exact  for  linear  elastic  behaivor.  Near  free  edges,  the 
boundary  conditions  of  the  LPT  solution  are  that  the  resultant  force  and  bend- 
ing moments  on  the  edges  are  zero.  The  LPT  model  permits  nonzero  stresses  on 
the  laminate  free  edges  while  the  true  boundary  conditions  are  zero  normal  and 
shear  stresses  on  these  edges. 

In  this  study,  the  tangential  stress  distributions  around  a circular 
hole  for  several  laminates  was  obtained  by  the  LPT  model  and  compared  with 
the  tangental  stress  distributions  from  the  three-dimensional  finite  element 
solutions  obtained  by  the  SAP  IV  program.  The  LPT  solutions  were  obtained  by 
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evaluating  the  effective  modulii  of  the  laminate  and  using  these  properties 
in  a two-dimensional  solution  for  an  orthotropic  plate  containing  a hole. 

These  properties  are  given  in  Table  A-l  of  Appendix  A.  The  stress  distribu- 
tions obtained  from  the  orthotropic  plate  analysis  are  resultant  stresses. 

The  stresses  in  each  ply  can  then  be  evaluated  from  the  LPT  distributions  by 
using  the  conditions  that  the  inplane  strains  are  the  same  in  all  plys.  Com- 
parisons of  the  tangential  stress  distributions  around  circular  holes  were 
made  with  results  obtained  from  the  three-dimensional  finite  element  stress 
analysis  and  results  from  the  LPT  model.  Twelve  laminates  and  the  three- 
dimensional  results  are  given  in  Tables  A-l  thrugh  A-8.  A summary  of  these 
results  and  the  LPT  results  are  shown  in  Tables  A-12  through  A-14.  Note 
that  the  LPT  solution  is  independent  of  stacking  sequence  and  stresses  do  not 
vary  through  the  thickness  of  the  ply.  This  is  not  the  case  for  the  three- 
dimensional  finite  element  results.  Thus,  a plane  had  to  be  selected  in  the 
three-dimensional  model  for  comparison  with  the  LPT  results.  The  midplane  of 
the  laminate  was  selected  as  the  location  for  obtaining  the  three-dimensional 
results.  The  reason  for  this  was  that  the  tangential  stresses  at  this  location 
appeared  to  take  on  the  largest  magnitudes.  A summary  of  comparisons  for  aQ 
given  in  Tables  A-12,  A-13,  and  A-14.  Note  that  the  stresses  for  the  (+8)  ply 

were  evaluated  from  (O-/+0/+0)  or  (9O„/+0/+0)  laminates  while  the  stresses 

4 — s 2 — s 

for  the  (07)  and  (90  ) plys  were  evaluated  from  (+0/+0/Oo)  and  (+0/+0/9CL) 
laminates.  The  greatest  variation  between  these  two  results  occurs  in  the 
(0y)  plys  of  the  (+0/+0/O  ) laminates.  The  differences  between  the  results 
for  the  (90^)  and  (+0)  plys  in  the  (+0/+0/9O2)g  laminates  is  about  the  same 
magnitude. 

Two  additional  laminates  were  considered  in  this  study  to  further 
illustrate  the .comparison  of  three-dimensional  finite  element  results  and  LPT 
solutions.  The  lay-up  is  denoted  by  (0,/+45/0)  and  (+45/CL/0)  . For  this 
laminate,  the  properties  used  for  the  (+45)  plys  and  the  0°  plys  ai  given 
in  Table  3.  Note  that  the  (+45)  plys  were  represented  by  as  a single  material 
with  no  extensional/bending  coupling.  Also,  the  value  for  of  .778  is  less 
than  what  would  be  obtained  by  rotating  the  properties  for  the  0°  ply  given 
in  Table  3.  Nonetheless,  these  properties  have  some  value  in  comparing  two- 
and  three-dimensional  solutions.  Figure  3 shows  a comparison  of  the  tangential 
stress  distributions  obtained  from  the  three-dimensional  and  two-dimensional 
analyses  for  the  (OjA+AS/O^  laminate.  These  differences  are  of  the  same  order 
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FIGURE  3.  STRESS  DISTRIBUTIONS  AROUND  A HOLE 

IN  A (0-/+45/0)  LAMINATE 
l — 8 


w 
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TABLE  3.  MATER IAL_ PROP ERT I ES  USED  IN  MODELING 

(0./+45/0)  AND  (+45/0o/0)  LAMINATES 
L — s — is 


Modulus 

Poisson's 

Shear 

Ratios 

Modulus 

(10  psi) 

£ 

(10  psi) 

0°  Plys  E2 


(+45)  Plys  E2 


30.0 

2.7 

2.7 

V12  ” 
V13  " 
V23  ’ 

.21 

.24 

.25 

G12  - 1.06 
G13  - 1.06 
G23  - 1.06 

3.77 

V12  “ 

.778 

G12  * 7*92 

3.77 

V13  “ 

.21 

G13  - 1.06 

2.7 

V23  " 

.21 

G23  - 1.06 

as  those  shown  in  Tables  A-12,  A-13,  and  A-14.  Figure  4 is  of  more  interest 
in  that  it  shows  how  the  tangential  stress  distribution  through  the  laminate 
thickness  and  compares  this  with  the  LPT  solution.  The  most  pronounced 
difference  appears  to  be  that  the  stress  "concentration"  for  the  0°  plys  is 
higher  in  the  three-dimensional  analysis  than  in  the  two-dimensional  analysis. 


5.  Effects  of  laminate  Thickness  to  Hole  Radius  Ratio 


In  the  case  of  a straight-free  edge,  the  only  factor  changing  the 
inplane  stresses  is  the  free  edge.  For  the  case  of  the  circular  hole,  the 
in  plane  stresses  are  alao  influenced  by  the  stress  concentration  due  to  the 
hole.  In  the  two-dimensional  solution  for  plate  containing  a circular  hole, 
the  region  of  influence  depends  on  the  hole  radius  and  the  material  properties. 
Also,  the  gradient  of  the  inplane  stress,  oQ,  is  very  steep  along  the  direction 
perpendicular  to  the  direction  of  loading.  Thus,  several  distinct  character- 
istics of  the  circular  free-edge  problem  do  not  appear  in  the  straight  free- 
edge  problem.  For  this  reason,  it  is  not  obvious  that  the  circular  free  edge 
behaves  like  the  straight  free  edge  even  for  small  laminate  thickness  to  hole 
radius  ratios.  Thus,  the  ratio  of  laminate  thickness  to  hole  radius  becomes 
an  important  parameter  to  study  in  obtaining  a basic  understanding  of  the 
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FIGURE  4.  DISTRIBUTION  OF  TANGENTIAL  STRESS  THROUGH 
LAMINATE  THICKNESS  AT  0 - 90° 


behavior  around  a circular  hole.  An  other  reason  for  examining  the  influence 
of  this  parameter  is  related  to  the  feasibility  of  obtaining  numerical  solutions. 
The  finite  elements  must  not  exceed  a certain  aspect  ratio  to  maintain  numerical 
stability  of  the  solution  technique.  Thus,  modeling  individual  boron  or  graphite 
laminas,  which  are  .005  to  .008  inch  thick,  could  require  more  finite  elements 
than  is  feasible.  However,  for  thick  laminates  where  the  lamina  thickness  to 
hole  radius  ratio  is  .5,  the  aspect  ratio  requirements  can  be  met  with  a feasible 
number  of  elements.  Thus,  if  the  solution  for  a thick  laminate  is  easier  to 
obtain  or  the  only  solution  available,  it  is  important  to  understand  how  the 
thick  plate  is  related  to  the  solutions  for  a thinner  laminate. 

To  help  understand  this  relationship,  a laminate  containing  a 1-inch 
hole  was  considered.  This  laminate  consists  of  two  materials  denoted  by  PLYA 
and  PLYB.  The  lay  up  for  this  laminate  is  (PLYB2/PLYA,./PLYB2) . Properties 
for  these  two  materials  are  given  in  Table  4.  The  hole  size  was  kept  constant 
and  the  thickness  of  the  laminate  was  allowed  to  vary  from  1.8  inches  to  .45 
inch. 


TABLE  4.  MATERIAL  PROPERTIES  FOR  THE  STUDY  OF  EFFECTS  OF 
LAMINATE  THICKNESS-TO-HOLE  RADIUS  RATIO 

(plyb2/plya5/plyb2) 


Extensional 

Poisson's 

Shear 

Ply 

Modulus 

Ratio 

Modulus 

Designation 

(psi) 

(psi) 

E1 

m 

30.0  x 106 

V12  “ 

.25 

G12 

- .65  x 106 

PLYA 

E2 

m 

2.7  x 106 

V13  “ 

.25 

G13 

- .65  x 106 

E3 

m 

2.7  x 106 

V23  " 

.25 

G23 

- .5  x 106 

E 

1 

m 

2.08  x 106 

V12  " 

.598 

G12 

- 2.38  x 106 

PLYB 

E2 

m 

2.08  x 106 

V13  " 

.136 

G13 

- 5.65  x 106 

E3 

m 

2.7  x 106 

v23  " 

.136 

G23 

- 5.65  x 106 
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Because  of  symmetry,  only  one-eighth  of  the  laminate  was  considered 
in  the  finite  element  representation.  This  is  one  quarter  of  the  top  half  of 
the  laminate.  Thirty-six  elements  were  used  in  this  model.  Three  elements 
through  the  thickness  and  twelve  elements,  as  shown  in  Figure  1,  were  used 
in  the  plane  of  the  laminate.  The  loading  conditions  were  presecribed  stresses 
on  the  boundaries  to  represent  the  laminated  plate  theory  solution  for  a uni- 
form average  stress  of  unity  in  the  1 direction.  The  applied  stresses  are 
given  in  Table  5. 

TABLE  5.  STRESS  BOUNDARY  CONDITIONS  FOR  THE 
(plyb2/plya5/plyb2)  LAMINATE 


Material 

o2  (psi) 

(psi) 

t12  <Ps1> 

PLYB 

.033 

.137 

0. 

PLY  A 

-.026 

1.690 

0. 

Four  cases  were  run  with  thickness  to  diameter  ratios  of  1.80,  1.35, 
.9,  and  .45.  The  tangential  stress  distribution  around  the  circular  hole  at 
the  laminate  midplane  is  shown  in  Figure  5 for  a ratio  of  0.9.  The  tangential 
stress  distribution  for  the  other  three  cases  did  not  vary  by  more  than  .98 
which  is  10  percent  of  the  maximum  value  at  9 * 90  shown  in  Figure  5.  The 
values  for  the  maximum  cr  decreased  with  decreasing  T/D  to  a value  of  8.8  for 
T/D  a equal  to  .45.  The  tangential  stress  distribution  at  the  midplane  was 
selected  because  the  largest  values  of  o occurred  there.  The  results  for 

U 

az  were  more  sensitive  to  the  changes  in  thickness.  Figure  6 shows  the  dis- 
tribution of  around  the  circular  hole  at  the  laminate  midplane  for  the  four 

ratios  of  thickness  to  hole  diameter,  T/D.  While  the  changes  in  o were  small 

z 

compared  to  the  maximum  value  of  the  tangential  stress,  the  changes  were  larger 

compared  to  the  magnitude  of  o . As  can  be  seen  from  Figure  6,  the  tensile 

z 

value  of  decreased  by  30  percent  as  the  ratio  of  T/D  changed  from  1.80  to 
.45.  For  this  cae,  the  value  of  ranges  from  -28  to  +16  percent  of  the 

average  applied  stress. 
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FIGURE  5.  STRESS  DISTRIBUTION  AROUND  A 

CIRCULAR  HOLE  AT  LAMINATE  MIDPLANE 
FOR  0°  CENTER  PLY 
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FIGURE  6.  DISTRIBUTION  OF  INTERLAMINAR  STRESS  Oz  AT 
LAMINATE  MIDPLANE  AROUND  A CIRCULAR  HOLE 
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SECTION  IV 


COMPARISON  WITH  EXPERIMENTAL  RESULTS 

In  this  section,  a comparison  of  experimental  data  and  results  from 
the  three-dimensional  finite  element  model  and  the  LPT  model  are  presented 
for  a (02/+45/0)s  laminate.  The  tangential  strain  distribution  around  the 
circular  hole  was  selected  for  the  comparison.  Data  were  obtained  from 
References  [12]  and  [21].  The  properties  in  the  material  direction  used  in 
the  analysis  are  given  in  Table  3.  In  the  three-dimensional  finite  element 
model,  the  (+  45°)  plys  were  treated  as  a single  material.  The  tangential 
strain  distribution,  en,  was  obtained  and  divided  by  the  far-field  strain  in 
the  direction  of  the  loading.  The  experimental  results  were  obtained  using 
the  measured  strain  concentration  factors  in  Reference  [12]  at  0°  and  90°. 
Data  for  intermediate  angles  were  based  on  the  birefringence  data  distribu- 
tions in  Reference  [12]  for  the  lowest  loading  level.  In  addition,  the 
tangential  strain  distribution  based  on  the  LPT  model  was  obtained.  The  com- 
parison of  predicted  and  measured  strain  distributions  is  shown  in  Figure  7. 
This  figure  shows  good  agreement  between  all  four  strain  distributions.  It 
also  shows  that  the  tangential  strain  distribution  in  the  (+45°)  plys  is  not 
the  same  as  that  in  the  0°  plys.  However,  the  average  is  close  to  the 
values  obtained  from  LPT. 
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FIGURE  7.  TANGENTIAL  STRAIN  DISTRIBUTIONS  AROUND  A 

CIRCULAR  HOLE  OBTAINED  FROM  THREE-DIMENSIONAL 
FINITE  ELEMENT  MODEL,  LPT  MODEL  AND  EXPERIMENTAL 
DATA  FOR  A (0J+45/0)  1AMINATE 


SECTION  V 


APPROXIMATE  ANALYSIS  FOR  THE  SIGN  OF  THE 
INTERLAMINAR  NORMAL  STRESS  AT  MIDPLANE 

While  the  three-dimensional  solutions  provide  a better  understanding 
of  the  interlaminar  stress  distribution  around  holes,  one  of  the  reasons  for 
conducting  this  study  is  to  obtain  solutions  which  can  serve  as  a basis  for 
developing  and  evaluating  simpler  stress  analysis  models.  In  this  section,  one 
such  model  is  presented.  This  model  is  very  preliminary  in  the  sinse  that  it 
is  applied  to  a limited  number  of  problems  and  based  on  certain  assumptions 
about  the  stress  behavior  of  the  laminate.  Nonetheless,  this  simple  model 
serves  as  a starting  point  for  attacking  a complex  problem.  The  limitations 
and  assumptions  of  the  models  are  described  in  the  following.  The  model  is 
then  applied  to  the  stress  loaded  laminate  problems  to  predict  the  sign  of  o 
at  selected  points  around  the  free  edge.  Good  quantitative  agreement  was  ob- 
tained between  the  three-dimensional  finite  element  results  and  the  results 
of  the  model. 

The  model  is  based  on  an  equilibrium  description  of  a wedge  shaped 
section  of  plate  shown  in  Figure  8.  Loading  conditions  are  stresses  applied 
to  the  boundaries  of  the  plate.  Unlike  the  model  of  Pagano  and  Pipes  [7],  this 
model  is  not  based  on  the  far  field  stresses  in  the  loaded  plate.  The  conceptual 
basis  for  this  model  is  that  the  Interlaminar  stresses  around  a circular  hole 
can  be  viewed  as  a reaction  to  a set  of  stresses  equal  to  but  having  the  op- 
posite sign  of  the  edge  stresses  from  the  Laminated  Plate  Theory  solution. 

This  concept  is  implemented  in  the  following  three  steps.  The  first  step  is  to 
obtain  a Laminated  Plate  Theory  (LPT)  solution  for  the  problem  of  laminated 
plate  containing  a circular  hole  with  stress  boundary  conditions.  The  LPT 

solution  gives  nonzero  stresses,  o (6)  and  t _(8),  for  each  ply  around  the 

r rt) 

hole.  The  second  step  is  to  impose  the  negative  of  these  edge  stresses  on  an 
otherwise  stress-free  body  and  determine  their  effect  on  the  interlaminar  stresses 
through  equilibrium  considerations.  The  third  and  final  step  is  to  superimpose 
these  two  solutions  to  obtain  an  approximate  solution  for  the  three-dimensional 
problem  LPT  stresses  around  the  circular  hole. 

Several  assumptions  concerning  the  equilibrium  of  the  section  shown 
in  Figure  8 are  worthwhile  to  pointing  out  before  writing  the  equilibrium  equa- 

21 
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tlons.  First,  only  two  equations  are  considered.  These  are  statements  that 

the  sum  of  the  forces  in  the  z direction  are  zero  and  the  sum  of  the  moments 

in  the  0 direction  are  zero.  The  shear  stress,  , is  not  included  and 

0z 

equilibrium  of  moments  in  the  r direction  is  not  considered  in  this  simple 

model.  The  assumed  form  for  o (r)  is  shown  in  Figure  8.  Due  to  the  sim- 

z 

plicity  of  this  model,  the  results  are  treated  in  a qualitative  way  here  and 

only  the  sign  of  o^  at  the  radius  of  the  hole  is  considered.  Since  is  the 

only  stress  contributing  to  force  equilibrium  in  the  z direction,  the  sign  of 

a at  the  edge  of  the  hole  can  be  found  from  the  sign  of  the  term  in  the  moment 
z 

equilibrium  equation  that  is  due  to  a r).  As  a result  of  the  preceding  as- 
sumptions, equilibrium  equations  for  the  force  in  the  z direction  and  the  moment 
in  the  0 direction  are  given  by 


0,+A0  R 

/ / o (r)  rdrd0  * 0 (1) 

' ' z 

0-0,^  r = a 

and 

0.+A0  R . 

/ / o (r)  (r -a)  rdrd0  - a (a)  — {h  + h,}aA0 

0-0,  r -a  z r z 1 z 

6j+A0  R 

♦ //  -5?  h2ilhi+h2ldrde  ■ 0 • <2) 

0=0^  r - a 

where  R is  selected  as  the  distance  from  the  edge  of  the  hole  where  the  free- 

edge  effects  are  negligible.  The  quantities  h,,  h2>  and  a are  shown  in  Figure  8. 

The  procedure  is  to  evaluate  the  second  and  third  terms  of  Equation 

(2)  using  LPT.  The  sign  of  o^  at  the  edge  of  the  hole  was  obtained  from  the 

sign  of  the  first  term  in  Equation  (2).  That  is,  if  the  first  term  in  Equation 

(2)  is  positive,  then  the  o is  negative  or  compressive  at  the  edge  of  the  hole. 
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Several  examples  are  considered  using  this  model.  Results  are  com- 
pared with  the  three-dimensional  finite  element  solutions  and  are  shown  in 
Table  6.  Results  for  two  locations  on  each  laminate  are  shown  in  this  table. 
These  locations  are  at  the  edge  of  the  hole  at  6 * 0,  and  6 - 90°.  Note  that 
there  is  agreement  between  this  model  and  three-dimensional  finite  element  re- 
sults for  6 ■ 90®.  However,  at  0 * 0°,  there  is  some  disagreement  for  the 
laminates  containing  (902>  plys.  There  could  be  several  reasons  for  this. 

One  possible  source  of  this  disagreement  is  that  the  assumptions  of  the  simple 

model  are  violated  in  these  cases  and  the  model  should  for  example  include  t . 

0Z 

While  this  model  is  simplistic  in  nature,  it  does  provide  a basis  that  can  be 

expanded  to  include  more  representative  stresses  and  equilibrium  equations. 

An  attractive  point  about  this  model  is  that  the  sign  of  a at  the  edge  of 

z 

the  hole  is  governed  by  the  stresses  from  the  LPT  solution  which  are  easily 
obtained. 


TABLE  6.  COMPARISON  OF  SIGNS  FOR  THE  INTERLAMINAR  STRESS  o 
OBTAINED  FROM  THE  THREE-DIMENSIONAL  ANALYSIS  AND  Z 
AN  APPROXIMATE  ANALYSIS.  T - Tension,  C - Compression 


Laminate 

6-0° 

0 

- 90° 

Approximate 

Analysis 

Three-Dimensional 

Analysis 

Approximate 

Analysis 

Three-Dimensional 

Analysis 

(+30/+30/02)s 

C 

-.203 

T 

.149 

(♦45/+45/0-) 

l s 

C 

-.345 

T 

.186 

(+60/+60/0.) 

l 8 

C 

-.376 

T 

.161 

(0./+30/+30) 

2 — 8 

T 

.068 

C 

-.081 

(0./+45/+45) 

T 

.181 

C 

-.062 

(0./+60/+60) 

2 ~ s 

T 

.233 

* 

C 

-.086 

(+30/+30/90o) 

C 

.051 

* 

T 

.255 

(+45/+45/90.J 
^ 2 s 

C 

.025 

T 

.659 

(+60/+60/902)a 

C 

-.142 

T 

.716 

(90./+30/+30) 

2 — 8 

T 

* 

-.114 

* 

C 

-.190 

(90./+45/+45) 

2 *—  8 

T 

-.135 

* 

C 

-.332 

(902/+60/+60)s 

T 

-.066 

C 

-.262 

* Denotes  disagreement  between  sign  of  o from  Approximate  Analysis  and  Three- 
Dimensional  Analysis. 
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SECTION  VI 


SUMMARY  AND  CONCLUSIONS 

Three-dimensional  finite  element  solutions  for  the  stress  distribu- 
tions around  a circular  hole  in  a laminated  plate  were  presented.  Problems 
were  selected  to  address  questions  concerning  the  effects  of  stacking  sequence, 
lay-up  angle,  and  a uniform  temperature  change  on  the  inplane  and  interlaminar 
stress  distributions.  A comparison  of  stresses  obtained  by  the  LPT  model  and 
the  three-dimensional  finite  element  model  was  made.  The  effects  of  the  ratio 

of  laminate  thickness  to  hole  diameter  on  the  interlaminar  stress  a was 

z 

examined.  Predicted  and  measured  strain  distributions  around  a circular  hole 
was  compared  and  a preliminary  investigation  of  an  approximate  model  to  predict 
the  sign  of  was  described. 

The  stress  analyses  of  the  laminates  were  obtained  using  the  SAP  IV 
three-dimensional  finite  element  stress  analysis  computer  program.  Conclusions 
are  thusly  based  upon  the  laminates  selected  and  the  results  of  the  SAP  IV 
finite  element  program. 

The  effect  of  stacking  sequence  on  the  interlaminar  stress  distribution 
around  a circular  hole  was  investigated  using  (0/90)  , (90/0)  , (O-/+6/+0)  , 

S 8 l — S 

(+0/+6/O-)  , (90 „/+6/+0)  , and  (+0/+0/9O„)  laminates.  Values  of  0 equal  to  30, 

— ZB  L — 8 — L 8 

45,  and  60  degrees  were  selected.  Loading  was  uniform  stress  in  the  0°  direction 
and  a uniform  temperature  change.  The  results  from  the  SAP  IV  three-dimensional 
finite  element  model  showed  that  changing  the  stacking  sequence  changed  the 
magnitude  and  sometimes  the  sign  of  the  midplane  interlaminar  stress  distribu- 
tion for  o z around  the  hole.  The  sign  of  o changed  for  some  stacking  sequences, 
but  no  general  rule  was  apparent  for  the  relationship  between  the  sign  of 
and  the  far  field  stresses.  The  (0/90)g  and  (90/0)g  comparisons  for  stress 

loading  showed  a change  in  sign  of  the  o stress  distribution  and  also  stress 

z 

change  in  magnitude.  For  example,  a changed  from  .25  for  the  (90/0)  laminate 

z s 

to  .03  for  the  (0/90)  laminate  at  6 - 90°.  Based  on  the  results  for  the 

s 

(O2/+0/+0)g  and  (90 ^1 ’+0 /+0) fi,  there  was  a trend  for  larger  values  of  a to  occur 
in  the  laminates  with  the  (Oj)  and  (9O2)  plys  in  the  center.  The  sign  of  the 
0^  distribution  for  the  (+60 /+60 / 90 2)  laminate  has  a sign  change  at  0 * 90’  for  a 
modification  of  the  stacking  sequence  which  places  the  ^Oj)  plya  on  the  out- 
side, However,  at  0 - 0,  the  sign  of  0 is  unchanged  for  a change  in  stacking 
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sequence.  The  effect  of  a uniform  temperature  loading  condition  coupled  with 
.)  change  in  stacking  sequence  was  that  the  (+30/+30/902)g  and  (902/+30/+30)g 
.iminates  showed  a change  in  the  sign  of  distributions  around  the  hole.  The 
other  two  laminates  had  o distributions  that  did  not  change  sign  around  the 
circular  hole.  However,  the  sign  and  magnitude  of  these  stress  distributions 
changed  with  stacking  sequence. 

A study  of  the  effects  of  lay-up  angle  in  the  (O2/+8/+0) g , (±6/+e/0^) 

(90^/+6/+6) g,  and  (+0/+0/9O2)g  laminates  for  values  of  0 equal  to  30,  45,  and 

60  degrees  was  presented.  Two  types  of  loading  conditions  were  considered;  an 

applied  stress  in  the  0°  direction  and  a uniform  temperature  change  to  represent 

the  fabrication  process.  The  results  of  the  uniform  stress  loading  applied  to 

the  (+0/+0/O-)  and  (+0/+0/9CL)  laminates  were  that  the  value  of  a on  the  edge 
2 s — 2 s z 

of  the  hole  at  an  angle  of  90°  to  the  loading  increased  as  the  lay-up  angle  in- 
creased from  30  to  60  degrees.  Becasue  of  the  constant  stress  loading  conditior 
and  the  reduction  in  compliance  of  the  laminates  for  increasing  values  0,  there 
is  a higher  applied  stress  in  the  0°  ply  for  increasing  0 as  shown  in  Table  A-2 
This  means  that  the  inplane  tangential  stress  in  the  (O2)  plys  and  the  (9O2) 
plys  at  0 * 90°  also  increases  as  0 increases.  The  maximum  tensile  value  of  o 

increased  for  increasing  values  of  8 for  the  (+0/+0/O-)  , (+0/+09O_)  , and 

— 2 s 2 s 

(O-/+0/+0)  . The  temperature  loading  condition  gave  a maximum  value  of  a for 
z s z 

the  (9O2/+0/+0)g  laminate  with  a lay-up  angle  of  0 = 45°. 

The  study  directed  at  the  effect  of  residual  stresses  due  to  fabri- 
cation of  the  laminate  was  based  on  a temperature  change  of  -200  F.  Values  of 

the  interlaminar  stress  0 from  the  tensile  range  of  800  to  1820  psi  resulted. 

z 

Inplane  tensile  stresses  were  as  high  as  9750  psi.  This  indicated  that  the 
residual  stresses  can  be  an  important  factor  in  predicting  the  response  of 
laminated  composites,  particularly  failure  modes  such  as  dedamination  or  ply 
splitting. 

A comparison  of  the  LPT  stress  distributions  and  those  obtained  from 
the  three-dimensional  finite  element  model  showed  that  there  was  not  close 
agreement  for  the  tangential  stresses  around  a circular  hole.  It  cannot  be 
concluded  that  all  laminates  would  give  similar  results.  The  closeness  of  the 
LPT  and  three-dimensional  results  appears  to  be  related  to  the  difference  in 
Poisson's  ratios  of  the  plys.  For  example,  results  were  obtained  (but  not 
shown  in  this  report)  for  Poisson's  ratios  of  .3  and  .5  for  adjacent  groups  of 
plys  indicated  close  agreement  between  LPT  and  three-dimensional  tangential 
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stresses  around  the  circular  hole.  The  comparisons  showing  the  larger  differences 
between  LPT  and  three-dimensional  results  occurred  for  Poisson's  ratios  of  .3 
and  1.2  for  adjacent  groups  of  plys. 

The  effect  of  laminate  thickness  was  considered  for  a particular 
laminate  with  stress  loading  conditions.  Ratios  of  laminate  thickness  to  hole 
diameter  in  the  range  of  .45  to  1.80  were  considered.  For  the  cases  considered, 
the  highest  value  of  tangential  stress  occurring  around  the  circular  hole  changed 
by  about  10  percent.  The  distribution  for  the  interlaminar  stress,  o^,  maintained 
the  same  shape,  but  decreased  in  magnitude  as  the  ratio  of  laminate  thickness  to 
hole  diameter  decreased.  In  general,  the  peak  tension  or  compression  stresses 
at  a ratio  of  laminate  thickness  to  hole  diameter  of  .45  were  about  2/3  of  the 
values  at  a ratio  of  1.8. 

A comparison  of  the  LPT  and  three-dimensional  finite  element  results 
with  experimental  data  was  done  for  one  laminate.  The  comparison  showed  good 
agreement  between  the  predicted  and  measured  tangential  strain  distributions 
around  a circular  hole.  This  provides  a degree  of  confidence  in  the  three- 
dlmenslonaal  model.  While  these  results  indicate,  that  the  inplane  behavior 
around  the  circular  hole  can  be  predicted  with  LPT  and  are  in  agreement  with 
the  three-dimensional  results,  it  is  pointed  out  that  the  Poisson's  ratios  for 
adjacent  plys  were  close  (.3  and  .78)  and  that  the  study  Is  too  short  to  be 
generalized  for  all  laminates.  For  example,  while  the  average  strains  around 
a circular  hole  appeared  close  to  the  LPT  solution,  the  strain  distribution 
varied  from  the  0°  plys  to  the  + 45°  plys  as  well  as  within  each  ply. 

The  results  of  this  study  show  that  the  interlaminar  stress  distri- 
butions around  a circular  hole  have  different  characteristics  than  thos  of  the 
straight  free  edge.  It  appears  that  the  stress  concentration  and  the  variation 
in  stress  distributions  around  the  circular  hole  contribute  to  the  complexity  of 
the  problem  of  relating  the  sign  of  to  the  field  stresses.  Thus,  equilibrium 
alone  may  not  be  sufficient  to  determine  the  sign  of  o^.  Consideration  of  other 
information  such  as  compatahility  may  be  needed.  For  example,  in  the  straight 
free  edge  problem.  References  [6]  and  [18]  showed  that  the  interlamlnate  stress 
oz  between  the  first  and  second  plys  of  a four-ply  system  could  not  always  be 
predicted  from  equilibrium.  This  stress  could  be  of  a different  sign  than 
that  at  the  midplane.  Similarly,  in  the  study  presented  here  for  circular  holes, 
the  sign  of  the  Interlamlnate  stress,  o^,  off  the  midplane,  was  not  always  the 
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same  as  that  at  the  midplane.  This  is  an  important  point  to  make  because  an 

understanding  of  the  signs  of  o ^ is  needed  if  one  is  to  consider  a relationship 

between  the  magnitude  of  o and  the  initiation  of  the  delamination  mode  of  failure. 

z 

One  important  result  of  a three-dimensional  study  of  the  stress  dis- 
tributions around  holes  in  laminated  plates  is  to  provide  information  and  results 
to  develop  and  gain  confidence  in  a simpler  model  that  can  be  used  routinely. 

As  a first  step  toward  accomplishing  this,  an  approximate  analysis  for  a wedge 
of  a laminate  was  used  to  estimate  the  sign  of  o ^ at  selected  points.  The  input 
to  this  model  Is  the  distribution  of  stresses  around  the  hole  obtained  from  the 
LPT  solution.  Comparisons  of  results  from  this  modt 1 and  the  three-dimensional 
finite  element  results  showed  agreement  in  19  out  of  the  24  cases  examined. 
Agreement  between  the  approximate  analysis  and  the  three-dimensional  finite 
element  results  occurred  at  0 “ 90  degrees  from  the  applied  load.  There  was 
agreement  in  all  cases  considered  at  this  location.  Disagreement  occurs  between 
the  approximate  analysis  and  the  three-dimensional  analysis  occurred  at  0 = 0° 
from  the  load  direction  for  5 out  of  12  cases.  If  the  far  field  stressed  would 
have  been  used  at  this  location,  disagreement  would  have  occurred  in  7 our  of 
the  12  cases.  While  the  approximate  analysis  presented  ignores  the  interlaminar 
shear  stresses  in  the  equilibrium  equations,  the  model  can  be  expended  to  in- 
clude these  stresses. 

In  concluding,  the  three-dimensional  finite  element  models  presented 
here  have  been  selected  to  address  several  important  questions  pertinent  to  the 
design  and  analysis  of  laminated  plates  containing  a hole.  It  is  emphasized 
that  only  a small  number  of  cases  were  considered  here  and  the  conclusions  are 
based  on  these  cases.  Within  this  context,  the  study  does  provide  insight  into 
the  stress  behavior  around  holes  in  laminated  plates  and  a means  for  constructing 
as  approximate  analysis  to  evaluate  the  sign  of  the  interlaminar  stress, 
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APPENDIX 

PROPERTIES  AND  RESULTS  FOR 
(+e/+e/o2)g  and  (+e/+e/902)s  lay-ups 


APPENDIX 


PROPERTIES  AND  RESULTS  FOR 

(+e/fe7o7j  and  (+0/Te/9O„)  lay-ups 

L 8 — L 8 


The  purpose  of  this  appendix  is  to  present  the  material  proper- 
ties, stress  boundary  conditions,  and  stress  distribution  for  the  analyses 
of  (+0/+0/O,)  , (O,/+0/+0)  , (+0/+0/9O)  .and  (90 ,/+0/+0)  . 

For  each  by-up  values  of  9 equal  to  30°,  45°  and  60°  were  con- 
sidered. Material  properties  and  stress  boundary  conditions  are  given  in 
Tables  A-l  and  A-2  respectively.  Results  for  both  mechanical  stress  and 
thermal  loading  conditions  are  presented  in  Tables  A-3  through  A-8  and 
Tables  A-9  through  A-ll,  respectively. 
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TABLE  A-l.  MATERIAL  PROPERTIES^  FOR  (+e/+0/02)g  LAMINATES 


Material  Property 

Ply8 

°2 

(+30/+30) 

(+45/+45) 

(+60/+60) 

E j (106psi) 

22.0 

8.54 

3.44 

2.03 

E2  (106psi) 

1.6 

2.03 

3.44 

8.54 

E3  (106psi) 

1.6 

1.65 

1.68 

1.65 

V12 

.25 

1.24 

.72 

.29 

V13 

.25 

-.035 

.08 

.18 

v23 

.25 

.18 

.08 

-.035 

o1?  (106psi) 

1.0 

4.54 

5.73 

4.54 

an  (106psl) 

1.0 

.88 

.78 

.70 

c23  (106psi) 

.64 

.70 

.78 

.88 

aL  (10_6/F) 

-1.5 

-2.99 

-.36 

6.04 

a2  (10'6/F) 

12.4 

6.04 

-.36 

-2.99 

a3  (10_6/F) 

12.4 

14.36 

15.3 

14.36 

(a)  Values  obtained 

from  Reference 

[19]. 
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TABLE  A-2 . STRESS  BOUNDARY  CONDITIONS  FOR  (+0/+0/O  ) 
AND  (+0/+O/9O2)s  LAMINATES  TO  PRODUCE  A^  S 
UNIT  STRESS  IN  THE  0°  OR  y DIRECTIONS 


Laminate 

(op 

or  (90°) 

Plys 

(+0/+6)  Plys 

a 

X 

0 

y 

(psi) 

o 

xy 

a 

X 

0 

y 

(psi) 

T 

xy 

(+30/+30/0,) 

-.094 

1.612 

0.0 

.047 

.693 

0.0 

(+45/+45/0o) 

- 1 s 

-.069 

2.243 

0.0 

.035 

.378 

0.0 

(+60/+60/02)s 

-.007 

2.531 

0.0 

.004 

.234 

0.0 

(+30/+30/90,) 

-.695 

.171 

0.0 

.348 

1.414 

0.0 

(+45/+45/90,) 

- 2 s 

-1.431 

.350 

0.0 

.716 

1.324 

0.0 

(+60/+60/90,) 

- 2 s 

-1.524 

1.702 

0.0 

.7620 

1.148 

0.0 
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TABLE  A-3.  STRESSES  AT  LAMINATE  MIDPLANE  FOR 

(+30/+30/02)s  AND  (+45/+A5/02)s  LAY-UPS 
WITH  UNIT  STRESS  IN~THE  Y-DIRECTION 


Parentheses  denote  stresses  from  adjacent  element. 


SIAidE  • 


(a)  Parentheses  denote  stresses  from  adjacent  element. 


TABLE  A-5.  STRESSES  AT  LAMINATE  MIDPLANE  FOR 

(+45/+45/902)8  AND  (+60/+60/902>g  LAY-UPS 
WITH  UNIT  STRESS  IN  THE  Y-DIRECTION 


(a)  Parentheses  denote  result  from  adjacent  element. 


TABLE  A-6.  STRESSES  AT  LAMINATE  MIDPLANE  FOR 

(02/+30/+30)s  AND  (02/+45/+45)s  LAY-UPS 


(a)  Parentheses  denote  result  from  adjacent  element. 


TABLE  A-7.  STRESSES  AT  LAMINATE  MIDPLANE  FOR 

(02/+60/+60)s  AND  (902/±30/+30)s  LAY-UPS 
WITH  UNIT  STRESS  IN  THE  Y-DIRECTTON 
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(a)  Parentheses  denote  result  from  adjacent  eleoent. 


TABLE  A-8.  STRESSES  AT  LAMINATE  MIDPLANE  FOR 

(902/+45/+45)s  AND  (902/+60/+60)s  LAY-UPS 
WITH  UNIT  STRESS  IN  THE  Y-DIRECTION 


(a)  Parentheses  denote  result  from  adjacent  element. 
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(a)  Parentheses  denote  result  from  adjacent  element. 


(a).  Parentheses  denote  result  from  adjacent  element. 


TABLE  A-ll.  STRESSES  AT  MIDPLANE  FOR  (+60/+60/902)s 
AND  (902/+60/+60)  LAY-UPS  WITH  THERMAL 
LOADINGS  OF  AT  - 200  F 


m 

m 

co 

cO 

On 

CO 

ON 

ON 

CM 

00 

0 

O 

r-^ 

r>» 

HT 

r*. 

r>. 

in 

CM 

m 

in 

O' 

cO 

O 

O 

NO 

CM 

m 

<r 

CM 

•n 

ON 

vO 

ON 

co 

m 

On 

m 

1 

H 

rH 

1 

1 

rH 

1 

1 

1 

• 

H 

ON 

00 

00 

CM 

00 

00 

r"» 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

v£> 

rH 

CO 

CO 

m 

<* 

O 

O 

NO 

sC 

O 

© 

s* 

O' 

vO 

00 

NO 

NO 

ro 

NO 

nO 

NO 

'T 

oo 

rH 

O' 

rH 

co 

H 

rH 

1 

1 

rH 

1 

1 

i 

CTJ 

r-N 

/— s 

r*v 

r“s 

vO 

° 

NO  rH 

rH  00 

CM 

O' 

rH  00 

r-  cm 

00  rH 

r-.  m 

no 

• • 

* • 

• • 

• 

• 

• • 

• • 

• • 

• • 

in 

P^  XT, 

nO 

oo 

m 

CO  CO 

ON  rH 

O rH 

no  m 

O 

co 

vO  rH 

no  ao 

ON 

vO 

00  CM 

m rH 

rH  CO 

00 

oo  oo 

m m 

On  On 

m f>*. 

co  ^ 

m 

m 

no 

p*.  on 

CO  CO 

O iO 

oo  oo 

•o  <r 

N— ' 

rH  »H 

rH 

rH 

rH  rH 

1 1 

1 1 

1 1 

« 

1 

» 1 

N-' 

v_^ 

N^ 

V 

S-x 

co 

rH 

m 

oo 

n. 

00 

CM 

NO 

o 

00 

CM 

CO 

rH 

00 

rH 

00 

00 

Ht 

to 

O' 

m 

o 

ON 

nO 

nO 

rO 

o 

vO 

rH 

ON 

vO 

ON 

00 

<5- 

*H 

rH 

1 

rH 

1 

rH 

1 

1 

rH 

CM 

/■n 

CC 

/■“N 

rs 

/-N 

r-s 

r-\ 

rv 

r\ 

/-v 

m rH 

CM  CM 

NO 

co 

m so 

no  m 

CM  O' 

rs. 

co  cm 

m m 

co  co 

co  o 

00 

ON 

m o 

rH  O 

m 

r-  vo 

O CO 

oo  oo 

vO  c^ 

vO 

O ^ 

00  o 

o 

rH 

NO  On 

CO  O' 

oo 

no 

co  co 

nO  On 

CO  vO 

CM 

O'  o 

CM  m 

m 

m 

o o 

cm  m 

in  in 

r- 

CM  rH 

CM  rH 

cm  cm 

rH 

rH 

rH  rH 

CO  CM 

CM  CM 

cn  cn 

N-r 

1 

1 

1 1 

N-*- 

i i 

1 i 

l i 

>-/ 

>— r 

w 

■w 

m 

<t 

CO 

ON 

ON 

NO 

m 

•H 

C30 

co 

NO 

o 

rH 

ON 

00 

>5- 

00 

rH 

o 

a 

m 

rH 

vO 

NO 

in 

NO 

o 

oo 

00 

co 

CM 

CM 

•O 

r* 

CM 

ON 

o 

*H 

NO 

1 

rH 

1 

NO 

CM 

1 

1 

CO 

1 

r*. 

rH 

H 

nO 

m 

ON 

Ht 

On 

m 

o 

co 

vO 

00 

CO 

CO 

O 

CM 

CM 

CM 

O 

O 

O 

CO 

nO 

co 

nO 

m 

CM 

NO 

in 

r- 

n* 

ON 

1 

rH 

1 

ON 

CM 

1 

1 

CM 

» 

w 

w 

0) 

X 

>N 

N 

8? 

«a J 

X 

N 

s? 

<T» 

u 

b 

b 

b 

H 

b 

b 

b 

b 

H 

D 

C/J 

» 

00 

«— s. 

CM 

© 

O 

Np 

O' 

1+ 

\ 

o 

O 

.* 

?, 

o 

CM 

2» 

O 

ON 

44 


I 


(a)  Parentheses  denote  result  from  adjacent  element. 


(a) 

TABLE  A- 12.  THREE-DIMENSIONAL  FINITE  ELEMENT  RESULTS  AND  LPT 
SOLUTIONS  FOR  TANGENTIAL  STRESS  DISTRIBUTIONS  FOR 
(02/+30/+30)8,  (+30/+30/02)s,  (0  /+45/+45)  AND 
(+45/+45/02)  LAMINATES 


e 

(+30/+30/02)b 
(02)  Plys 

(02/+3O/+30)s 
(+30)  Plys 

Three-Dimensional 

Analysis 

LPT 

Three-Dimensional 

Analysis 

LPT 

0° 

-.3884 

-.2835 

-.3630 

-.4637 

13.3° 

-.3258 

-.2242 

-.4746 

-.4646 

26.6° 

-.1250 

-.1026 

-.2050 

-.3892 

45.0° 

-.5612 

.0815 

.6258 

.3384 

63.4° 

1.5167 

1.7981 

1.7874 

2.1158 

76.7° 

4.8991 

4.5609 

2.2896 

2.6173 

90.0° 

7.2951 

5.8391 

2.4048 

2.5129 

(+45/+45. 0o) 

— Z 8 

(02/+45/+45)s 

(°2> 

Plys 

(+45)  Plys 

Three-Dimensional 

Three-Dimensional 

e 

Analysis 

LPT 

Analysis 

LPT 

0° 

-.3846 

-.2317 

-.5922 

-.8749 

13.3° 

-.3073 

-.1631 

-.7034 

-.8251 

26.6° 

-.0959 

-.0551 

-.1773 

-.5195 

45.0° 

-.7164 

.2739 

.8909 

.9065 

63.4° 

1.8534 

3.0262 

1.2629 

2.0203 

76.7° 

6.1819 

5.6950 

1.197 

1 . 5066 

90.0° 

9.3488 

6.7456 

.9826 

1.1379 

(a)  Evaluated  at  laminate  midplane. 
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TABLE  A- 13.  THREE-DIMENSIONAL  FINITE  ELEMENT  RESULTS^  AND  LPT 

SOLUTIONS  FOR  TANGENTIAL  STRESS  DISTRIBUTIONS  FOR 

(+60/+60/02)a,  (02/+60/+60)s,  (+30/+30/902)s,  AND 

(90-/+30/+30)  LAMINATES 
z — s 


e° 

(+60/+60/02)8 
(02)  Plys 

(0-/+60/+60) 

Z — s 

(+60)  Plys 

Th  ree-Dimens ional 
Analysis 

LPT 

Three-Dimensional 

Analysis 

LPT 

0° 

-.1715 

-1.1902 

-1.4142 

13.3° 

-.1226 

- .9620 

-1.1211 

26.6° 

-.0056 

-.0384 

- .1285 

- .2778 

45.0° 

-.6820 

.6484 

.8530 

1.1756 

63.4 

2.1189 

3.5333 

.8207 

1.5303 

76.7 

6.8667 

6.3164 

.6745 

1.0242 

vO 

O 

o 

10.3716 

7.5932 

.5980 

.7033 

(+30/+30/902)8 

(90^30 /+30)s 

(902) 

Plys 

(+30)  Plys 

Three-Dlmens ional 

Three-Dimensional 

e 

Analysis 

LPT 

Analysis 

LPT 

0° 

-3.2425 

-2.5310 

-.1551 

-.2344 

13.3° 

-2.1903 

-1.7857 

-.3373 

-.2895 

26.6° 

- .5072 

- .3183 

.1587 

-.1378 

45.0° 

.5650 

.6486 

1.5615 

1.1756 

63.4° 

.4734 

.4266 

3.0350 

3.0837 

76.7° 

.5845 

.4338 

3.7326 

3.9655 

0 

O 

.6312 

.5142 

3.6397 

4.2426 

(a)  Evaluated  at  laminate  midplane. 


TABLE  A- 14.  THREE-DIMENSIONAL  FINITE  ELEMENT  RESULTS (a)  AND  LPT 
SOLUTIONS  FOR  TANGENTIAL  STRESS  DISTRIBUTIONS  FOR 
(+45/+45/902)  , (902/+45/+45)  , (+60/+60/90  ) AND 
(90»/+60/+60)S  LAMINATES  8 

Z “ 8 


(+45/+45/90,) 

Z 8 

(902/+45/+45) 

s 

(902) 

Plys 

(+45) 

Plys 

Three-Dimensional 

Three-Dimensional 

6 

Analysis 

LPT 

Analysis 

LPT 

0° 

-5.0800 

-3.3959 

- .2646 

- .5728 

13.3" 

-3.4731 

-2.3695 

- .4984 

- .6268 

26.6° 

- .9730 

- .2040 

.7941 

- .1362 

45.0° 

.9052 

.7001 

3.0358 

2.3171 

63.4° 

.8656 

.3590 

3.0363 

3.3851 

76.7° 

1.0930 

.6357 

2.4902 

3.2162 

90° 

1.2216 

.8158 

2.2123 

3.0808 

(+60/+60/9CL) 

Is 

(902(+60/+60) 

8 

(902) 

Plys 

(+60) 

Plys 

Three-Dimensional 

Three-Dimensional 

e 

Analysis 

LPT 

Analysis 

LPT 

0* 

-5.7912 

-3.9941 

-1.1889 

- 1.7187 

13.3° 

-3.9022 

-2.3192 

- .7533 

- .13309 

26.6° 

-1.0050 

.4136 

1.1154 

.4867 

45.0° 

1.3842 

.7232 

3.1828 

3.0020 

63.4° 

1.4759 

. 7859 

2.5596 

2.9797 

76.7° 

1.8010 

1.2314 

1.9582 

2.5517 

90° 

1.9957 

1.4457 

2.0243 

2.3647 

(a)  Evaluated  at  laminate  midplane. 
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